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1^ 

' Abstract We derive a sequence of measures whose corresponding Jacobi matrices have 

special properties and a general mapping of an open quantum system onto ID semi infinite 
chains with only nearest neighbour interactions. Then we proceed to use the sequence of 
measures and the properties of the Jacobi matrices to derive an expression for the spec- 
tral density describing the open quantum system when an increasing number of degrees 
^ , of freedom in the environment have been embedded into the system. Finally, we derive 

' convergence theorems for these residual spectral densities. 

^3 ; 



1 Introduction 

, 1.1 Background 

, Real quantum mechanical systems are never found in complete isolation, but invariably 

Xf^ • coupled to a macroscopically large number of "environmental" degrees of freedom, such 

as those provided by electromagnetic field modes, density fluctuations of the surrounding 
media (phonons) or ensembles of other quantum systems, like electronic or nuclear spins 
[TK2113]. The fact that the environment is totally or partially inaccessible to experimental 
probing in these so called " open quantum systems" leads to the appearance of an effectively 
irreversible dynamics for the quantum system's observables, and mediate the fundamental 
.J processes of energy relaxation, phase decoherence and, possibly, the thermalization of the 

rS I subsystem. 

' Accurate numerical or analytical description of general open quantum systems dynamics 

appears, prima facie, to be extremely difficult due to the large (often infinite) number of 
bath and system variables which need to be accounted for. When the environmental de- 
grees of freedom are modeled as a bath of harmonic oscillators exact path integral solutions 
can be available but are rarely of practical use [TUIITII^. Hence assumptions such as weak 
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system-environment coupling and vanishing correlation times of the environment, i.e. the 
Born-Markov approximation, are often invoked to obtain compact and efficiently solvable 
equations. These approaches suffer the drawback that their accuracy is hard to certify and 
that they become simply incorrect in many important situations. Indeed, our increasing 
ability to observe and control quantum systems on ever shorter time and length scales 
is constantly revealing new roles of noise and quantum coherence in important biological 
and chemical processes [HOEUTUH] and requires an accurate but efficient description of the 
system-environment interaction that go well beyond the Born-Markov approximation [8j|9] 
in order to understand the interaction of intrinsic quantum dynamics and environmental 
noise. In many biological, chemical and solid-state systems, deviations from strict Marko- 
vianity, which can be explicitly quantified I1 HI121[T5] . are significant and methods beyond 
standard perturbative expansions are required for their efficient description. A number of 
techniques have been developed to operate in this regime. Those include polaron approaches 
[13], the quasi-adiabatic path- integral (QUAPI) method T5|, the hierarchical equation of 
motion approach |16j and extensions of the quantum state diffusion description to non- 
Markovian regimes [T7] . 

Here we will focus on the exploration of the mathematics of an exact, analytical map- 
ping of the standard model of a quantum system interacting with a continuum of harmonic 
oscillators to an equivalent model in which the system couples to one end of a chain nearest- 
neighbour coupled harmonic oscillators, as illustrated in (a) of Fig. ([1]) and (b) of Fig. ([Ij. 
This mapping has permitted the formulation of an efficient algorithm for the description of 
the system-environment coupling for arbitrary spectral densities of the environment fluc- 
tuations [18lll9j . This new mapping which originally intended just as a practical means of 
implementing t-DMRG which would avoid approximate determination of the chain represen- 
tation using purely numerical, and often unstable, transforms j20j, it was quickly realised 
that the scope of the mapping is much broader. Indeed, the mapping itself provides an 
extremely intuitive and powerful way of analysing universal properties of open quantum 
systems, that is, independent of the numerical method used to simulate the dynamics [T5] . 
This conclusion was implied also in [21 , where the authors also developed a, in principle, 
rather different chain representation of a harmonic environment using an iterative propa- 
gator technique [31] . 

Both of these theories establish chain representations as a novel and direct way of looking 
at how energy and correlations propagate into the environment in real time. In the chain 
picture the interactions cause excitations to propagate away from the system, allowing a 
natural, causal understanding of Markovian and non-Markovian dissipation in terms of the 
properties of the chain's couplings and frequencies. An intuitive account of the physics of 
chain representations, non-Markovian dynamics and irreversibility for the method in [191 
[18' is given in |22| , and an interpretation of chain parameters in terms of time-correlations 
for the method in [2] can be found in |23) . 

1.2 This communication 

The main goals of this paper can be split into two groups. Firstly, we aim at developing 
a general framework for mapping the environment of an open quantum system onto semi 
infinite ID chain representations with nearest neighbour interactions where the system only 
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couples to the first element in the chain. In these chain representations, there is a natural 
and systematic way to " embed" degrees of freedom of the environment into the system (by 
" embed" , we mean to redefine what we call system and environment by including some of 
the environmental degrees of freedom in the system) (c) of figure ([T|) . One can make a non- 
Markovian system-environment interaction more Markovian by embedding some degrees 
of freedom of the environment into the system, a technique already employed in certain 
situations in quantum optics |24j and that recently has been demonstrated in all generality 
[21j . What remains unclear however, is to quantify how efficient such procedures can be as 
well as determining the best way of performing the embedding. We will show that those 
issues can be efficiently addressed with the formalism presented in this manuscript. 

This general formalism allows the comparison between different chain mappings. Thus 
the second aim is to develop a method for understanding how Markovian such embed- 
dings are by finding explicit analytical formulas for the spectral densities of the embeddings 
corresponding to the different chain mappings. Furthermore, we also derive universal con- 
vergence theorems for the spectral densities corresponding to the embedded systems and 
give rigourous conditions for when these limiting cases are achieved. This paves the way for 
a system interacting with a complex environment to be recast by moving the boundary of 
the system and environment, so that the non-trivial parts of the environment are embedded 
in the new effective "system" and the homogeneous chain represents the new, and much 
simpler, "environment" - See (c) in figure ^ for a pictorial illustration. The advantage of 
this is that the residual part of the environment might be simple enough for some of the 
approximations mentioned in section 11.11 to be applied, enabling us to integrate out these 
modes and dramatically reduce the number of sites of the chain that have to be accounted 
for explicitly. In order to achieve this, we have to first develop new mathematical tools and 
theorems regarding secondary measures and Jacobi operators, greatly extending and devel- 
oping the application of orthogonal polynomials that was used in the original chain mapping 
of 18J. These results might be useful also in other areas of mathematics and mathematical 
physics which are related to the theory of Jacobi operators such as the Toda lattice |25) . 

We show that this generalised chain mapping reduces to two known results mentioned in 
section [Ol under special conditions. For the first known result [T51IT5] . the general method 
developed in this paper gives analytical and non-iterative expressions for the spectral den- 
sities corresponding to the embeddings. For the other special case [H], we derive calculable 
conditions for when the spectral densities corresponding to the embedding converge, -an 
aspect not addressed in |21| . As seen in the examples, we apply this technique to derive 
exact solutions for the family of Spin-Boson models which will allow us to illustrate how 
the different embedding methods are related. In addition, the method developed here is 
also valid when the spectral density of the system-enviroment interaction has a gap in its 
support region. This is of practical interest as there are open quantum systems (such as 
photonic crystals [26) ) that naturally exhibit such a spectral density and hence can only be 
mapped onto a chain using the method presented in this paper. 

The contents of each section is as follows: Section [5] is concerned with deriving the 
necessary mathematical tools for the application to open quantum systems in the subse- 
quent sections. We start by introducing some elementary results in the field of orthogonal 
polynomials in section 12.11 (this section also helps to introduce notation) . In section 12.21 we 
focus on deriving a formula which makes explicit a sequence of secondary measures solely 
in terms of the initial measure and its orthogonal polynomials. We point out that although 
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authors such as Gautschi introduce the concept and definition of secondary measures, here 
we provide for the first time an analytical closed expression for them in terms of the initial 
measure and its orthogonal polynomials. This result will turn out to be a vital ingredient 
in the development of the subsequent chapters. Moreover, Gautschi states that the general 
solution we have found is unknown f57^ pl6-17. In section [^31 we study the properties of the 
3-term recursion coefficients of the orthogonal polynomials generated from their measures. 
We then define beta normalised measures which have more general and useful properties. 
The main result of this section is our new theorem regarding the Jacobi matrix. This new 
theorem will be used extensively in section [31 Section [3] is concerned with chain mappings 
for open quantum systems and embeddings of environmental degrees of freedom into the 
system. In section [3. II we develop a general framework for mapping open quantum systems 
which are linearly coupled to an environment onto a representation where the environment 
is a semi-infinite chain with nearest neighbour couplings and define two special cases of 
particular interest, - the particle and phonon mappings. In section 13.21 we investigate the 
relation between this work and recent work by [21 . In order to do this we make extensive 
use of the relations developed in section 2. We show how their mapping is a special case 
of the work presented here and find analytical non iterative solutions to quantities such as 
the sequence of partial spectral densities. Section 13.31 is dedicated to deriving the formulas 
for the sequence of partial spectral densities for the particle case mapping. In section 13.41 
we develop convergence theorems for the sequence of partial spectral densities. We show 
rigorously that the sequence converges under certain conditions and give the universal func- 
tions the spectral densities converge to for the particle and phonon cases. The conditions for 
which the sequences converge are stated in terms of the initial spectral density. In section 
[4] we give explicit analytic examples for the family of spectral densities of the Spin-Boson 
model for the particle and phonon mapping cases. 

The main new results of this article are theorems (HH), <^5^, (HH), (El), (0, (HH), (ESJ and 
corollaries dH) and Q. 

2 Secondary measures 

2.1 Introduction to notation and basic tools 

Definition 1 Let us consider a measure dfi{x) = jl{x)dx with real support interval I heg- 
gining at a and ending at b which has finite moments: 



with > over I such that Co{d^) > 0. 

We note that throughout this article, unless stated otherwise lower and upper values a and 
b can be finite or infinite. 

Definition 2 Let P denote the space of real polynomials. Then, for any u{x) and v{x) G P 
we will define an inner product as 




(1) 




(2) 
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Definition 3 We call {P„(d/i; x)}^o ^^''^^ orthonormal polynomials with respect 

to measure where each polynomial P„ is of degree n, if they satisfy 

{Pn{dtl),Pm(dll))f:,= 6nra n, m = 0, 1, 2, . . . . (3) 

Similarly, 

Definition 4 We call {7r„((i/x; x)}^q the set of real monic polynomials with respect to 
measure d^, where each polynomial 7r„ is of degree n if they satisfy 

Hn{dn]x) ^ Pn{dn\x)/an n 0, 1, 2, . . . , (4) 

where a„ = an{dfi) is the leading coefficient of Pn{dfi;x). 

Theorem 1 For any measure dfi{x), there always exists a set of real orthonormal polyno- 
mials and real monic polynomials. 

Proof See [17]. 

Theorem 2 The monic polynomials 7r„ satisfy the three term recurrence relation 

7r„+i((i/x;x) = (x - a„)7r„(d^;a;) - /3„7r„_i(d^;a;) n = 0,l,2,..., (5) 
■K-i{dn\x) := 0, (6) 

where 

f, N {x^ln{d^l),^Tn{d^l))f, 
a„ = a„(rf/i) = - — — — , n = 0,l,2,..., (7) 

f3n ^ Pnidfi) ^ T— n = l,2,3, .... (8) 

{TTn-l(dfl),TTn-l(dfJ,))fi 

Proof See [17]. 

Definition 5 We will define ^^{dii) by 

Po{dfi) = {TTo{d^),TTo{dfi))fi. (9) 

Corollary 1 f^oid^) ~ Co{dfi). 

Proof We note that from definition dU, 7ro(d/z;x) — 1 for all measures d^{x). Hence, 

l3o{dfi) = {no{dn),no{d^i))f^ = rf^(.T) = Co(d/i). (10) 

J a 

a 

Theorem 3 When the measure dji has hounded support, the a„ (d/i) and /3„ [d^) coefficients 
are hounded hy 

a<anidfi)<b n = 0,1,2,..., (11) 

< I3n{dfi) <max{a'^,b'^) n = 0,1,2,.... (12) 
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Proof See [37]. 

Theorem 4 The orthonormal polynomials Pn{d^; x) satisfy the three term recurrence rela- 
tion 

tnPn+l{dtl\x) = [x - Sn)Pn{dlJ.]x) -tn^lPn^l{d^l\x) 71 = 0,1,2,..., (13) 

P-i{dfi;x) := 0, Po{dfi;x) = l/y/MM, (14) 

where 

Sn = s„{dfi) = an{dfi) n = 0,1,2,..., (15) 



tn = tnidfJ.) = \/ I3n+i{dn) n = 0,1,2,.... (16) 

Proof See [17]. 

Definition 6 We will call Qnid/i; x) the secondary polynomial^ associated with polynomial 
Pn{dfi;x) defined by 

n M ^ /''' Pn{dn]t) - Pn{dfi;x) 

Qn[dfi;x)^ / dti{t), n = 0,1,2,... (17) 

Lemma 1 The polynomials Qn{d^; x), n = 1, 2, 3, . . . are real polynomials of degree n — I 
and Qq — 0. 

Proof Follows from writing P„((i/i; a;) in the form Pn{d^; x) — X]g=o ^i?^'' : using the identity 
fi -x'i = {t- x) YZZl'^ tPx'i-^-P, g = 1, 2, 3, ... to cancel the denominator in Eq. ^ 
and noting that Co(dfi) > 0. □ 

Definition 7 The Stieltjes Transformation of the measure dfi{x) = fi{x)dx is defined bu\28\ 
where z e C — [a. fe] . 

Theorem 5 If a measure dp{x) — p{x)dx has Stieltjes transformation given by 

S,iz) = z - C,{dp) - (19) 

with Co^dfi) — 1, then the secondary polynomials {(3n(d/i; x)}'^^i form an orthogonal family 
for the induced inner product of dp{x). 

Proof See [53]: or [3^ for a direct proof. 



^ also known as polynomial of the second kind. 
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2.2 Derivation oj the sequence oj secondary normalised measures 

Definition 8 For two measures dp{x) and d^(x) satisfying Eq. (|19p . we call dp{x) the 
secondary measure associated with d/j^x). 

Definition 9 We call the sequence of measures dfiQ, d/ii, (i/i2, • ■ • generated from a measure 
dfiQ by 

Sp^^^{z)^z-Ci{dp„)--\- n = 0,l,2,..., (20) 

dpn{x) — flnix)dx n — 0,1,2, . . . , (21) 

l^n{x) = y^^ n= 1,2,3,..., (22) 

the sequence of normalised secondary measures, where Co(d/io) = 1- 

Corollary 2 All measures in a sequence of normalised secondary measures have their zeroth 
moment equal to unity. 

Proof By taking zeroth moment of both sides of Eq. we find 

Co(dMn) = §^ = l n = l,2,3,.... (23) 

We also have that Co((i/io) = 1 by definition □ 
Lemma 2 

n-l 

Cn{dpm+l) = Cn+2{dpm) " Ci(dflm)Cn+l{dprn) ^ ^ C s{dpm+i)Cn- s(dfJ.m) (24) 

s=Q 

n,m — 0,1,2, .. . . 

Proof For simplifity we will prove Eq. (I24p for m = as the generalisation is trivial. By 
Taylor expanding Sp,g{z) and Sp-^{z) in x = 1/z we find 

^MoW = f:%§^ asz^oo, (25) 
^p.W = E^^ (26) 

From Eq. (IT^ we have 

Sp,{z)SM = {z- Clidpo)) S,A^) ~ 1. (27) 
Hence substituting Eq. (US]) and (gH) into Eq. d^T]) we find 

C„(dA*o)C„,(rfpi)x"+"+2 = ( - - Ci(dMo) ) E Cs{dfio)x'+^ - 1. (28) 



n.ni—0 ^ ^ s— 
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By comparing terms of the same power in x and taking into account Coidfio) = 1 we deduce 
Cn{dfio)Cm-n{dpi) = Cm+2{d^a) " Ci (d^o)C*,„+i (d^io) m = 0, 1,2, . . . . (29) 
By a change of variable in Eq. ([2^ we finally arrive at Eq. (IMl) . □ 

Lemma 3 A sequence of normalised secondary measures d^o, dfii, dfi2,---, dfXn, can he 
written as a continued fraction of the form 

Soiz) = ^ 

"0 

z - C'lfi 

«i 

Z - Cii : 



z - Ci, 



dn-1 

Z — Oi,„_i — 



Z — Ci,n — dnSn+l{z) 

(30) 

where we have introduced the shorthand notation Sn{z) '.— S'p^(z), Cn,s ■— Cn(d^s)i n,s — 
0, 1, 2, . . . ; dn := C2^„ - „, n = 1, 2, 3, . . and da := 1. 

Proof By evaluating Eq. (1^^ for n = 0, and taking into account the above definition of dn 
we see that d„ = C{dpn), n = 1, 2, 3, . . .. Using our new notation, Eq. (PO)) for the sequence 
of measures reads 

5„+i(^) = -^[z-Ci,„--^], n = 0,l,2,.... (31) 

Solving this for S'„(z) followed by repeated substitution gives us Eq. (|5(7| . □ 
Theorem 6 The following relations hold for the continued fraction Eq. pOp 

Q / N Un(-dnSn+l{z)) + , , 

So{z}^ — , , n = 0,1,2,..., (32 

Vn{-dnSn+l{z)) + U„+l 

with relations 

Un+l = (z - Ci^n)Un - dn-lUn-l, U„+l = (z - Ci,„ )u„ - d„_iU„_i , (33) 

and starting values 

-uo = 0, ui = 1, Wo = 1, vi=z-Co,i- (34) 

Proof These are elementary results from the theory of continued fractions (e.g. see section 
4, connection with continued fractions |31jV 



Lemma 4 



An+l -.^ Un+lVn - Vn+lUn = dodid2 ■ . -dn-i, n= 1,2,3,..., (35) 

Al = 1. (36) 
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Proof Using Eq. (j33p to substitute for Un+i and v„+i into Eq. psp . we find tlie relation 
An+i = dn-i^n- Using Eq. (p4)) to verify Eq. (1351) for the starting values, Eq. ((35)) follows 
by induction. □ 

Definition 10 A Pade Approximant for a function g of type q/p in the neighbourhood of 
is a rational fraction 

Fi.) ^ 2M, (37, 

with degree of Q < q, degree of P < p and g{z) — of order 0(2^+^+^) in the neighbour- 
hood of 0. 

For more details, see |31) . 

Theorem 7 F„(z) = H2±i(£) ^5 a Pade Approximant for So{z) of type n/{n + 1), n = 
0,1,2,.... 

Proof Using theorem we can write Sa{z) ~ as 

Through lemma ^ we see that An+i is independent of z. By Taylor expanding Sn+i{z) 
defined in definition ([7]) about x = 1/z, and remembering that Cold fin) = 1 7t,= 1,2,3..., 
we find using Eq. (p8)) that 



Sn+i{z) ^- +0{^) asz^oo n = 0,1,2,.... (39) 

z 

By induction, we see that u„+i and m„+i, n = 0, 1, 2, . . . given by Eq. ([55]) and (IMl) . are 
degree n + 1 and n polynomials in z respectively, both with leading coefficients equal to 
unity. Hence we conclude that 

w„+i(z) A„+id„ dQdi...dn , . 

Soiz) - ^^^^^^^ = ^2n+3 ^ ^2n+3 aS Z ^ OO 71 = 0,1,2,.... (40) 

Thus by definition (|T0| we conclude the proof. □ 

Lemma 5 Un{z) = XnQnidfiQ] z) and Vn{z) — A„P„(o?/io; z) n = 0,1,2,..., with A„ = 
l/a„ iwere a„ is defined in definition 

Proof From section 5.3.: moment problems and orthogonal polynomials (p213-220) of [5T| . 
and theorem ([7]) it follows that 

Mn+l(^) Qn+i{dfio; z) 



Vn+liz) Pn+l{dHQ]z) 

By observing the starting values, we also have that 



n = 0,l,2,.... (41) 



(42) 



voiz) Po{dfiQ;z)' 
Hence 

Un{z) ^ XnQnidfio; z) and Vniz) = XnPnidHo'yZ) n = 0,l,2,.... (43) 
Given that u„ and w„ have leading coefficients equal to 1, we must have A„ — l/a„ 
n=0,l,2,... . □ 
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Theorem 8 = a?i/a^+i = /3„+i(d/^o) n = 0, 1, 2, . . . . 
Proof Proceeding in the same way as in page (18) of [27j, we have 
Qn{dfio;z) 7„ 

^°(^)' P„(d^o;.) n = 0,l,2,.... (44) 

where 

7n = — / a: Pn(dMo;a;)dA^o(a^) = ,n . , ^ „ . , n = 0,1,2,.... (45) 

Noting that Po{dfio; x) / ao = 7ro(d/io; x) = 1 and that Co{dfio) — 1, Eq. pS)) tells us 70 = 1. 
Comparing Eq. (gi]) with Eq. (gD]) and (gT]), we deduce that 

, , , _ flo {PnWo),P^{dno))-^ . „ o 

< {Poidno),Po{dfio))f, 

By induction it follows 

^ On {Pn+l{d^io),Pn+lid^J.o))f, ^^-^38 (47) 

^n+l {Pn{dflo),Pn{dfJ,o))fL 

Due to definition Q, we see that {Pn{d^o) i Pn{dfio)) p. = 1 n = 0, 1, 2, . . . , hence 

n = 1,2,3,.... (48) 



"n+l 



From definition we see that Eq. (|47)) can be written in the form 

, _ {Trn+l{dfio),TTn+l{dfJ.o))p „ _ r, o ^Q^ 

(7r„((i^o),7r„(rfA*o))/i 
Hence, from definition ([5]) we conclude, 

d« = /3n+i(rfAio) n = 1,2,3,.... (50) 

For n = 1, Eq. gives us 

(7ri(d/io),7ri(rf/io))p ag 

"0 = 7 — — ^ — 7^ — — - Pi(,a/ioj. (olj 

{no{dfio),TTo[d^io))fi a\ 

□ 

Theorem 9 Measures d/i(x) = fl{x)dx can be calculated from their Stieltjes transform by 



K^) = -J— 1™ [SJx ~ ie) - Sn{x + ie)] . (52) 
Zm £-5-0+ 



Proof This result is known as the Stieltjes- Perron inversion formula. See [3 2) . 
Definition 11 We call Lp{d^]x) the reducer of measure dfj,{x). It is given by 

(p{d^; x) = lim [Sp,{x — ie) + Sp_{x + ie)] . (53) 
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The reducer allows us to write an explicit expression for the secondary measure associated 
with as follows. 

Theorem 10 For a measure dfi{x) with secondary measure dp{x), we have 



p{x) 



i^^^P^+TT^fi^x) 



(54) 



Proof See 0^. 

Definition 12 We define the functions Zn{x)eC n — 0,1,2, . . . as 

Zn{x) = ^'^^'^ +iTTp.nix) 71 = 0,1,2,..., 

where (pn{x) := ip{dfin',x). 

Lemma 6 The following recursion relation hold for Zn^i{x) 



71 = 0,1,2,. 



(55) 



(56) 



Proof We start by finding a relation between the Stieltjes transformation of a measure 
dfi{x), and the reducer of its associated secondary measure dp{x): By definition, we have 



(p{dp;x)— lim [Sp{x — ie) + Sp{x + ie)] . 

e— >0+ 



Using Eq. (fTOj) . we find 



ip{dp; x) = 2[x — Ci{dp)\ — lim 



- ie) + Sp,{x + ie) 



>o+ Sfi,{x — ie)Sp{x + it) 
Now using theorem © and definition (fTT|) . wc find that 

Sfi,{x ~~ ie) + Sfj,{x + ie) _ <(5(cJm; x) 



lim 



0+ Sp{x-ie)S-^{x + it) ;en^^.^2^2(^) 



Hence from Eq. ([55]) we get 

(p(rfp; a:) = 2 [x - Ci (d/i)] - -^t^^.,) 



ip{dp;x) 



t:'^p'^{x) 



(57) 



(58) 



(59) 



(60) 



Using the definition of a sequence of normalised secondary measures, definition ^ and the 
definition of dn in Lemma from Eq. (l60t we find 



<<5n+i(a:^) 



1 



2[a;-Ci,„] 



v« (a;) 



(2;) 



71 = 0,1,2, 



(61) 



4 ' ""^MriW. 

Similarly, we can also write theorem (jlOp for our sequence of normalised secondary measures 
using definition ^ and d„. We find 



/i„+i(a;) = 

U.77 



tlnix) 



, 71 = 0,1,2,.. 



(62) 



4 ' ^Vr>(2;). 

If we write Eq. ([55]) for Z„+i(a;) and then substitute Eq. (|5T|) and Eq. (p^ into the RHS 
we arrive at Eq. (j56l) . □ 
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Theorem 11 The following relations hold between Zo{x) and Zn+i{x) 

Zo[z) = — , --^ , n = 0,1,2,..., (63) 

were Un{z) = XnQn{dfJ-o; z) and Vn{z) = A„F„(d^o; z) n = 0, 1, 2, . . . , with A„ = l/a„, a„ 
defined in definition 

Proof By comparing Eq. (|56p with eq. (PT|) . we note that Zn{x) satisfies the same recursion 
relation as Sn{x). Hence theorem ([5^ readily applies to Eq. ([55| if we exchange Sn+i{x) 
with Zn+i{x) and 5*0(2:) with 2'o(a;). The relations between m„, i;„ and (5„, P„ are proven 
in lemma ([5]). □ 

We are now ready to state our first main theorem: 

Theorem 12 A sequence of normalised secondary measures starting from dfia: d/io, dfii, 
d^2,- ■ ■ , dfirn,- ■ ■ can be generated from the first measure in the sequence d/io by the formula 

a M - I ^ f64) 



n= 1,2,3,... , 
where the tn coefhcients are defined in definition Q. 
Proof After solving Eq. (|63p for Zn{x), we find 

ry f \ o,n-i Zo{x)Pn{dfio;x) ~ Qn{dfio;x) _i o q faK\ 

^n[X — r7 / \ D 71 ^ 7^ 71 T n — i,Z,6,.... [00 

o„d„_i Zo(x)Pn-i(dfio\x) ~ Qn-i(dfio;x) 
From theorem ([8]) we see that 



an-l/On = Kny/dn-1 n= 1,2,3,..., (66) 

where k„ is the sign of a„_i/a„. After using this relation to simplify the a„_i/a„(i„_i 
coefficient in Eq. (p5|) and substituting for Zo{x) using definition ([T^. we take real and 
imaginary parts to achieve 

. , s Kn /2o(a;) [F„_i((i/io;a;)(5„(d/io;a;) - P„(d/xo;x)(3„_i(d//o;a;)] .„„s 
l^n{x) = ^=7 ^ -2 (67) 

Vrf"-i (^p„_,(d^o;^)Mi^a^_Q„_^(d^o;2:)) +^2^§(x)p2_,(d^o;x) 

n = 1,2,3,... . Using the identities from lemmas (jj]) and ([S]), we note that 

Pn{d^io; x)Qn+iid^j,o; x) - Pn+i{d^iQ; x)Qn{d^io; x) = dodi . . . d„_i/A„A„+i (68) 

n = 1, 2, 3, . . . . Using theorem ([8|) and lemma ([5|) to write the RHS in terms of the a„'s, 
we find 

dodi . . .dn-i/XnKi+i =^ alan+i/on ^ al/Kn+i\/d^ n = 0,l,2,..., (69) 
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where in the last hne we have used Eq. (1551) . We now note that by definition Po(d/io; x) — 
flQ. We also have by definition ([3]), that {PQ{d^o',x),Po{d^Oix))fio = 1. Hence 

1 = {Po{dfio;x),Po{dfiQ:x))fi„ = alCo{dfio) = Oq. (70) 

Hence from Eq. ([TO)), we find that we can write Eq. (|S7)l as 

Mn(a^) = -J ^ -2 (71) 

''"-1 (p„_i(d^o;:r)£%El_Q„_i(d^o;:r)j +n^f4{x)P?^_,{dfio;x) 

n = 1,2,3,.... Finally, with the relations from theorems (|5]) and ^ we find dn-i = 
tl^Mf^o) n = 1,2,3,.... □ 



2.3 Derivation of the Jacobi matrix theorem 

Lemma 7 The an{dfi) and finidfi) coefficients defined in theorem ^ are invariant under 
a change of scale of the measure dfi{x), while the change in ^^{dfi) scales linearly: 

aniCdn) = an{dfj,) n — 0,1,2,... (72) 
MCd^i) = Mdfi) n= 1,2,3,... (73) 
MCdfi) = CMdfi) (74) 

where C > 0. 

Proof First we will show that 7r„(C(i^; x) — 7r„(d^; x) n = 0, 1, 2, 

From definition (|4|), we have that 

1 = (P„(d/i),P„(rfAi))p n = 0,l,2,.... (75) 

By multiplying and dividing by C we find 

1 = ( ^ : n = 0,l,2,.... (76) 

From definition we conclude 

P„(CdM;x) = ^^ii^ n = 0,l,2,.... (77) 

By observing definition ^ and Eq. (l77|) . we conclude 

QniCdfi) = anidfi)/VC n = 0,l,2,.... (78) 
Hence, from Eq. dTT]), dZH]) and definition (gj it follows 

TTn{Cdfi;x) = Pn{Cdfi;x)/aniCdfi) = P„((i^; x)/a„(d/i) = 7r„(d^;a;) (79) 
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n = 0, 1, 2, . . . . Now we can see how the Q!„ and /3„ coefficients change: 
Using definition ([2]) and Eq. (l79l) we have 

a„(Oc(/ij = - — , , , — (,»Uj 

(7r„ (Cd/X ),TTn{Cdfi))c-fj, 
^ {x7Tn{dlJ,),TTn{dlJ,))p. 

{Trn{dfi),nn{dfi))fi 
= n = 0,l,2,.... (82) 

Similarly, we find 

/3„(CdM) = /3„(d/i), n = 1,2,3,.... (83) 

In the case of /3o we have 

MCdfi) - {MCd^l),noiCd^l))^^ = CMdfi). (84) 

□ 

Theorem 13 If dp{x) is the secondary measure associated with d^[x), then 

a„+i(d^) = Q!„(rfp) rj = 0,1,2,... (85) 
/3„+i(dAi) -/3«(rfp) 71-1,2,3,... (86) 

Proof See [57], theorem 1.36 (page 16). 

Lemma 8 A sequence of normalised secondary measures dfj.o{x), d^i[x), d/i2(a;), . . satisfy 

am+n{d^J.c,) ^ ctn{d^im) n, m = 0, 1, 2, . . . , (87) 
/3„+„(d/xo) = /3„(d^m) n = 1, 2,3, . . . , m = 0, 1,2, . . . . (88) 

Proof Using lemma ([7]) and definition ([9]) we find 

anidpm) ^ anidpjn/Coidpm)) = an{d^im) n = 0, 1,2, ... TO 1,2,3 ... , (89) 
Pnidpm) ^ l3nidpm/Co{dp„i)) ^ I3nidfl,n) n = 1 , 2, 3, . . . TO = 1 , 2, 3 . . . . (90) 

Hence taking into account theorem (fT5)l wc find 

a„+i(d^p) = Q;„(d^p+i) = 0, 1, 2, . . . , (91) 

Pn+iidpp) = lin{dpp+i) 71 = 1, 2,3, . . . , p = 0, 1, 2, . . . . (92) 

We will now proceed to prove Eq. ([57]) by construction: by evaluating Eq. (pij) . for (n,p) 
at (s — 1, to), (s — 2, to + 1), (s — 3, TO + 2),. . . ,(0, to + s — 1) for any to > 0, s > 1 we have 
the following sequence of equations 

as{dp^-,) = as-i[dp,n+i) (93) 
as~i{dpm+i) = as-2{dp,n+2) (94) 
as-2{dpm+2) = as-z{dp,yi+z) (95) 



ai[dpra+s-i) ^ ao{dprn+s)- (96) 
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Thus by repeated substitution we arrive at 

Q;s(d/im) = ao(dAfm+s), m = 0,1,2,... ,s = 1,2,3,... . (97) 

We note that this equation is also vahd for s = 0. If we make the change of variable s = p+n, 
m = in Eq. (|97| followed by relabeling indices, we find 



am+s{d^J.o) ^ ao{dnm+s) m, s = 0,1,2, ... . (98) 

Hence from Eq. ^ and dH]), we arrive at dST]). Similarly, we prove Eq. □ 

Definition 13 We call the sequence of measures dvQ^x), di'i{x), dv2{x), dv3{x), ... beta 
normalised measures, where di>Q{x) defines the sequence 

-n{x) = — , (99) 

(Pn-i(rft/o;a:) '^ 2°"^ ^ Qn-i{dvo\x)\ + 7T^iy^{x)P^_j^{diyo; x) 

71 = 1, 2, 3, ... , and dvn{x) — Vn{x)dx n = 0, 1, 2, . . . . 

Lemma 9 For every sequence of beta normalised measures {dvn}^"!^ , there always exists 
a sequence of normalised secondary measures {d/i„}"^g° such that 

Vn{x) ^ Pn{dvQ)fLn{x) 71 = 0,1,2,.... (100) 

Proof For ti = in Eq. (llOOp and taking into account corollary ([Ij, we have 

vii(x) = Co{dvo)p.o{x). (101) 

The only additional constraint on dfiQ as compared with any measure dfi, is that Co((i/io) = 
1. By calculating the zeroth moment of both sides, we see that relation Eq. poip satisfies 
this additional constraint. Now we will proceed by finding an expression for the sequence of 
normalised measures generated from d/iQ in terms of dvQ. For this we need to find how the 
quantities Pn{d^o', x), Qnid^o; x), /3„+i(d/io), ti = 0, 1, 2, . . . and ip{d^io; x) can be written in 
terms of Pn^dvo; x), Qn{dvQ; x), (3n+i{dvQ), ii — 0,1,2, . . . and ip{dvQ; x) respectively. Using 
relation Eq. (jlOlD and definition Q, we find 



Pn{dti^;x) = ^C^{dvo)Pn{dvo;x) 71 = 0,1,2,.... (102) 
Using this relation, and definitions ^ and (ITT|) . we find 

Qn{diyo;x) 

Qy,[diio;x) = —=== 71 = 0,1,2,..., (103) 
VCo(rft'o) 

if{dvQ;x) 

ip{dfio;x) ^ (104) 
Co[dvo) 

Lemma ([7]) tells us 

/3„(d/io) =/3„(d7^o) 71=1,2,3,.... (105) 
Hence using relation Eq. (fT6)) . we find 

tl_Ml^o) = tl-iidi^o) = (3n{diyo) 1,2,3,.... (106) 
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Now substituting Eq. (fTU^ . (fTU5)) . (fTUi)) . and (|TM1) into Eq. 

1 Po(a;) 



(107) 



n = 1, 2, 3, . . . . Hence by observing definition p^ . we find Eq. (|100p for n = 1, 2, 3, . . . . 
□ 

Theorem 14 ^ sequence of beta normalised measures dvQ{x), dvi{x), dv2{x), . . ., satisfy 
ayn+n{dvo) ^ a-nidvyn) n, TO = 0, 1, 2, . . . , (108) 

Prn+nidvo) ^ PnidVm) n, TO = 0, 1, 2, . . . . (109) 

Proof From lemma ([9]) we see that dun and di^n are related by a constant, hence using 
lemmas Q and ([8]) we find 

Q;„i+„(di'o) = a„(c?t'm) = 0, 1, 2, . . . , (110) 

l3rn+n{diyo) ^ l3nidVrn) n = 1 , 2, 3, . . . , TO = 0, 1 , 2, . . . . (Ill) 

For (3a{dvm) we find 

;3o(dl^m) = /3o(/3m(rft'o)rfMm) = /5m(c?l'o)C'o(d/im) TO = 0,1,2,..., (112) 

where we have used Eq. (|100p followed by Eq. ([7^ and then corollary ([1]). But Co(fi^m) = 
1 TO = 0, 1, 2, . . . , by definition. □ 

Remark 1 As we will see in section SI for a wide range of dv^; am+nidvo) , /3m+„(di'o) and 
dvm can be determined analytically. Hence Eq. (llOSp and Eq. (|109D can also be used to find 
analytical solutions to a wide range of integrals. 

Definition 14 We will call the infinite tridiagonal matrix of a measure dfj,(x) 

01 



Jid^^) 



aoidfi) ^/Pi{d^i) 



ViMdfi) a2{dfi) ./Psidfi) 







(113) 



the Jacobi matrix. 5*66 |27j for more details. 
Definition 15 We will call the matrix 



Jnidfi) = 



an{dn) y^/3n+l{dfl) 
^/ I3n+l{d^i) an+i{d^i) l3n+2{dlJ.) 



yjl3n+2{dlj) a„+2(d/i) y/ /3n+3i.dfJ.) 







(114) 



n = 1, 2, 3, . . . , the nth associated Jacobi matrix of the Jacobi matrix J^{dfi). 
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We are now ready to state our second main theorem: 

Theorem 15 For Jacobi matrices for which its corresponding measure defines a sequence of 
normalised secondary measures, there exist an infinite sequence of associated Jacobi matrices 
corresponding to the sequence of normalised secondary measures. These matrices are formed 
by crossing-out the first row and column of the previous Jacobi matrix in the sequence: 

Jn{d^Ji^)=J{diJin) n=l,2,3,... (115) 

Proof By equating the matrix elements in Eq. (|115p . we find Eq. ((87)) and (|88)) . Hence 
lemma ® implies Eq. (fTTS)) . □ 

Corollary 3 Theorem (1151) is also valid for any sequence of measures which are proportional 
to a sequence of normalised secondary measures such as the beta normalised measures. 

Proof Given that the Jacobi matrix does not contain the /3o coefficient, the result follows 
easily from lemma ([7]). □ 



3 Chain mappings of open quantum systems and Markovian embeddings 

3.1 Generalised mapping 

An open quantum system can be represented by a system plus bath (also known as envi- 
ronment) model introduced by Caldeira and Leggett [23] ■ The Hamiltonian for this model 
is 

H = Hs + He + H,^t, (116) 

where Hs is the Hamiltonian which governs the system dynamics, He is the Hamiltonian 
of the environment and Hint describes the interaction between the two. In the case where 
the environment has a continuous set of degrees of freedom and the interaction between the 
system and environment is linear and does not depend on the state of the system, we can 
write the Hamiltonian as 

H = Hs + He + H„,t = Hs + dxg{x)ala^ + / dxAh{x){al + a^). (117) 

where g{x),h{x) G M and have real domain [xmin,Xmax\- A is a arbitrary hermitian op- 
erator acting on the system part of the Hilbert space. The set of operators {aa;,aj,}, x S 
[xmimXmax] are crcatiou and annihilation operators and hence satisfy [a2,,a|^] = 5{x — y). 
Physically g{x) represents the dispersion relation of the environment and h{x) determines 
the system-environment coupling strength. Together they determine the spectral density in 
the following way 

Definition 16 We call the function J{x) the spectral density, 

J(c.)^./.2[g-i(.)]^^^, (118) 

where g^^{g{x)) ~ g{g^^{x)) — x, Dom[J] e [0,uj„iax], and J{uj) > its domain. For cases 
where uJmax — oo, we adapt the definition of the domain in an obvious way. 
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Remark 2 As we will see later in the proof to theorem (ITSl) (more specifically in Eq. 11451) , 
we use the J{x) function in the definition of a measure M{x) in such a way that AI{x) = 
if J{x) = 0. Hence if we want to form a valid measure according to definition ([1]), we need 
J(x) > over its domain. For cases where this does not hold and J{x) is zero at a point(s) 
in the interval {Q^uj„iax), we can define a set of spectral densities by redefining their domain 
such that they contain only non positive values at the boundaries of their domain. Also see 
remark (p]). 

The influence of the environment on the reduced system's dynamics is uniquely determined 
by J{x) when the initial state of the environment is gaussian [ID] . Consequently, Eq. (|118p 
tells us that a specific choice of g{x) does not determine uniquely the system's dynamics 
when the initial state on the environment is gaussian. This freedom of choice of g{x) will 
be vital in the proof of theorem which is our third main theorem. We should also 

note that theorem (jl6p is also valid for any initial environment state if the g{x) and h{x) 
functions found in Eq. (|142p and Eq. (|144p respectively, are equal to the g{x) and h{x) in 
the initial Hamiltonian, i.e. Eq. (jll9p . 

Theorem 16 A system linearly coupled with a reservoir characterized by a spectral density 
J{x) is unitarily equivalent to semi-infinite chains with only nearest-neighbors interactions, 
where the system only couples to the first site in the chain. In other words, starting from 
an initial Hamiltonian, 

H = Hs + He + H,nt = Hs + dxg{x)ala.^ + / dxAh{x){al + a^), (119) 

there exists a unitary operator {?^n(a^)}5JLo such that the countably infinite set of new oper- 
ators 



bi= dxUn{x)al n = 0,1,2,..., (120) 

satisfy the corresponding commutation relations 

[bn,bl]^Snm n, m = 0, 1, 2, . . . , (121) 
with transformed Hamiltonian 

H ^Hs + He^, + H,nt,q , (122) 

OO 

He., = J2 {Ein{q){e^'^'bibi + e-2'^^6„&„) + E2n{q)bibn (123) 

n=0 

+EM{^'''^'^nbi+i + e-^'^-bnbn+i) (124) 

+£;4„(g)(6tfe„+i + 6„6l+i)}, (125) 

H,nt,q = E5n{q)A{e'^'bl + e-'^^bo). (126) 

where Ein(q), E2niq), E^n{q), E^niq), E^niq), A2 and q G [0, 1] is a free parameter of 
the mapping which determines the particular version. For a pictorial representation of this 
theorem, see (a) and (b) of figure ([T]). 
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Proof The proof is by construction. 

Let us start by defining the most general type of transformation of the creation and an- 
nihilation operators aj, and Ox in terms of another set, Cx and cj,, which preserves the 
commutation relations. We can do this via the so called Bogoliubov transformation^^ 

ax = e'^i^^) cosh(r(x))c^ + e*''^^^) sinh(r(a;))4, (127) 
4 = e-'^i^^) cosh(r(x))4 + e"^'^'''^ sinh(r(a;))c^, (128) 

where 6*1 (x), 02ix), r{x) G M and [a^jjaj^] — [cx,cl^] — 6{x — y)^liwe now parameterise r(x) 
by introducing a new function ^(x) through r{x) ~ In^(x) where ^(x) G [0, +oo) using Eq. 
(jl27p and (jl28p we find after renormalising 

dxg{x)aiax = dx^^({^\x) - l){e^^<'^^^ c^x^l + e^^^^^^cc.) 

+2{e{x) + l)clcx) (129) 



where AO{x) — 02 (x) — Oi{x). We note that as far as Eq. (jl29p is concerned, we can set 
9i{x) = —62{x) without losing any generality. With this definition of 6*1 (x), the system- 
environment term simplifies to 

A dxh{x){ax + al) = A dx/i(x)e(x)(e'^^(^)4 + e-*''^(^)c:,). (130) 



For appropriate choice of the function ^(x), we can define a measure 

dX{x) = M{x)dx (131) 
M(x) = h''{x)e{x). (132) 

This measure allows us to define a set of monic polynomials {7r„((iA; x)}^q, which we can 
use to define the set of functions {[/„(rfA; x)}5^q through the relation 



U [dX-x)^ ^ 7r„(dA;x)/t(x)g(x) ^^^^^ 

(7r„(dA),7r„(dA))M \/ (7r„(dA), 7r„(dA))M 
We can now define the set of creation and annihilation operators of the chain 

4 = / dx?7„(dA;x)4, (134) 



with inverse relation 



&„ = / dxUn{dX]x)cx, (135) 



t =^[/„(dA;x)6t, (136) 

oo 

. = ^f/„(dA;x)6„. (137) 



ri=0 



^ We can replace the cosh and sm/i functions with cos and sin functions for fcrmions. 



20 



M.P. Woods et al. 



Let us now assume that the phase factors are constant]^ ^2(2;) 
stitution of Eq. (|T36| and Eq. (fT37)) into Eq. ([130)) we obtain 



Ao^constant. After sub- 



A 



'bo). 



(138) 



So we note that for all functions ^(x) that result in a valid measure Eq. (|13ip . one can 
achieve a coupling between system and reservoir which only interacts with the first element 
in the chain. 

Now we will examine carefully what type of chain can be generated via this generalised 
mapping. 

Using the orthogonality conditions of the orthogonal polynomials and Eq. (jl36p and 
(jl37p . we can transform terms of the form CaxCbxdx into terms of the form J2'^=o W„banbbn 

where the WnS are constants and the sub indices a,b — denote that the operator is 
an annihilation operator and a,b — 1 denote that they are creation operators. Also, we 
can transform terms of the form J^^^"^ xCaxCbxdx into X]^o ^nbanbbn + Wiba(n+i)bbn + 
W2nbanbb{n+i) by usiug the three term recurrence relations Eq. ([5]) to eliminate the x. We 
can also map terms of the form J'j^'""^" x'^ CaxCbxdx fc = 2, 3, 4... using the three term recur- 
rence relations k times, but this would result in every chain site coupling to its fcth nearest 
neighbourslfl Given this resoning, we speculate that there is only one trial solution to Eq. 
(|129p using real polynomials with A9{x) — 292{x) — 2A2 which will map it onto a chain of 
nearest neighbour interactions. This is: 



2e{x) 



i^\x)+l)=C2+g2X. 



Solving these gives us 



ax) 



(C2 + 2ci) + (32 + 2gi)x 



(c2 - 2ci) + (52 - 2gi)x 



1/4 



gix) = y/[{c2 - 2ci) + (32 - 2.9i)a;][(c2 + 2ci) + (52 + 2gi)x] 



(139) 
(140) 

(141) 
(142) 



We will now work out the corresponding measure. Using the formula for the differentiation 
of an inverse function, 



dg'Hx) 



and Eq. (fTTS)) . we find 



dx 



h\x) 



(x) = 1/ 

Jigix)) 



dg{x) 



dx 

dg{x) 



dx 



ig-\x)) 



(143) 
(144) 



^ If we were using complex polynomials, we could leave it as a function of x and hence we would 
have to impose a less restrictive constraint. See footnote ((6]). 

* We could also transform terms of the form J^™'"" f{x)caxCbxdx using orthogonal polynomials in 
f{x), but this will result in the same chain mapping if we choose f{x) = x and using the standard 
orthogonal polynomials in x, so we will stick with the simple x case only. 
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so using Eq. (jl3ip . we have that the measure M{x) is 



M{x) = ft'(x)C'(a;) 



dgjx) 
dx 



{x)e{x) 



(145) 



Let us define a particular relation between the independent variables ci , C2 , gi , 52 , namely 



2.91 
2ci 



9.92, 

-9C2, 



where 9 G [0, 1]13 We now find using Eq. and (ITi5|) 



(1 - 9)c2 + (1 + 9).92a; 



-,1/4 



(1 + 9)c2 + (1 - q)92X 



9(x) = + 9)c2 + (1 - 9)322:] [(1 - 9)c2 + (1 + 9)522:], 

M(.)^j(,(.))gi (;+/;--+;;--^;)-'^-- . 

TT (1 +9)C2 + (1 - 9)522; 

For convenience, we will denote the quantities in Eq. (jl48l) . (jl49l) and (I150p by 



e(52x) 
5'(52:r) 
g2M^g2x) 



9{x), 
M{x), 



(146) 
(147) 



(148) 
(149) 
(150) 



(151) 
(152) 
(153) 



where ^^(2;), g'^{x) and M''{x) are independent of 52, to remind us of their q and 92 
dependency. We will also denote 



dA« = g2M'i{g2x)dx, 
Pn{q) = MdX")- 



(154) 
(155) 
(156) 



First we note from definition ([T5| that the support of (x) is determined by the domain of 
J(x). As we have chosen the domain of J{x) to be [0, oJmax], we need to choose the support 
of g2M'^{g2x) such that the J{x) is integrated over its full domain. By observing Eq. p50p . 
we see that we need to choose Xmin such that g''{xmin) — 0. There are two such values 
which satisfy this requirement, 



Xmin± = -(c2/ff2)(l ± 9)/(l T 9)- 



(157) 



We will take Xmin^ as it does not contain any singularities in q. Furthermore, we will set 
C2 = 9/4 as this will allow us to compare the extremal solutions of the mapping easily 
as we will see later. We note that for 52 > the measure has no singularities for well 
defined spectral densities J{x). Also ^"^(2:) and g'^ix) are positive functions in their domain 
[xmin- , Xmax] for all Xmax > Xmin^ ■ Furthermore, the measure dX' satisfies the requirements 



* In principle there are other choices for these constants which will result in other types of 
coupling relations between the chain elements. We have chosen this particular case as it will result 
in a particularly interesting solution. 
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of definition ([IJ and hence it is a valid measure. We can also find a value for Xmax- By 
observing Eq. (I150p we see that we need Xmax to satisfy g'^{g2Xmax) — '^max- There is only 
one solution to this which satisfies Xmax > Xmin- which is 

x„^ax ^(q{q^ + l)- 2v/g4 - 4c.^^Jg2 _ 1)) / (4^^ {q^ - 1)) . (158) 

For now on we will write Xmin- = Xmin^il) and Xmax = Xmax{q) to remind us of the q 
dependency. The transformation is now well defined with no degrees of freedom other than 
q and 52- 

Now let us perform the mapping. After substituting Eq. p36p and (|137[) . into Eq. (|129p 
and using the orthogonality conditions, we fin^fl 

/ dxg2g'^{g2x)ala^ 

Jx„,,„_(q) 

+ g2y/K^ {qie'^'^'bibi^i + e-2'^^6„6„+i) + (fet 6„+i + 5„6l+,)) }. 

□ 

Remark 3 Note that in the generalised mapping, theorem (jl6p , the spectral density can 
represent continuous modes and/or discrete modes. The discrete modes are represented by 
dirac delta distributions and the continuous modes by continuous functions. In the case of 
N discrete modes only, the generalised mapping will map the system onto a chain with N 
sites. Mathematically, the reason for this is because the inner product of the measure ^ 
in this case can only be defined for functions living in the space spanned by a set of N 
orthogonal polynomials of finite degree (See discrete measure, page 4 following Theorem 
1.8 of [27j). Physically, this is because the number of degrees of freedom corresponding to 
the Hamiltonian after and before the mapping have to be the same. Also see remark ([6]) 
regarding other features of spectral densities. 

In addition, it is also worth noting that if the system interacts linearly with more than 
one environment through different system operators A for each environment, then one can 
easily generalise the results of the generalised mapping (theorem (jl6p ) such that we can 
map the Hamiltonian onto a Hamiltonian where the system interacts with more than one 
chain (i.e. one chain for each environment). 

Lemma 10 Without loss of generality, we can set g2 = 1 in theorem (jl6p . 



We note that in order to perform the mapping, we have had to set the phase factors of the Bo- 
goliubov transformation 9i{x) and ^2(2;) equal to constants. From (jlSOp we see that this restriction 
means that they play the simple role of a change of phase of the new variables b„ and ftj,. If one 
were to work with complex polynomials, one could have probably imposed less stringent constraints 
and then the phase factors would play a more interesting role such as a gauge transformation. 



(159) 

(160) 
(161) 
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Proof First we will denote dX'^{x) = dX'^{g2,x) to make explicit their 52 dependency. From 
definition we have that the induced norm of the monic polynomials {7r„((iA''(52); a;)}^o 
is generated by the inner product 

{TTn{dX'>{g2);x),TTr,{dX''{g2);x))dxH92) n-0,1,2,.... (162) 
By making the change of variable y — x/g2, we can write this as 

{TTn{dX'^{g2);x),TTn{dX'^{g2);x))d\i(g2) = (tTk (dA« (1) ; g2y) , 7r„ (dA'' (1) ; ^2^)) dA? (1) (163) 

n = 0, 1, 2, . . . . However we note that this change of variable means that the polynomials 
TTnidX'^ (1); g2v) n = 0, 1, 2 in the RHS of Eq. (|163l) are monic polynomials in g2y and hence 
are not monic polynomials for measure dX'^{l,y). To make them monic polynomials in y, 
we must make their leading coefficient unity. Hence by dividing both sides of Eq. (|163p by 
gl" we get the relation 

/ > u\if-i\ w (7r„ (dA*? (ff 2 ) ; a:) , 7r„ (dA'? (ff 2 ) ; a:)) dAng2 ) 
(7r„(dA''(l);y),7r„(dA«(l);?;))rfA<!(i) = 2^ — (164) 

n = 0, 1, 2, . . . . Similarly, we find 

/ M\in\ \ ^^^q|^\ ^\ (a;7r„(dA'?(g2); x), 7r„(dA'?(ff2); a:))dA^(g2) 
(y7r„(dA'?(l);y),7r„(dA''(l);2/))rfA9(i) = 2Wri (^^^) 

92 

n = 0,l,2,... . Hence using Eq. pM)) and together with Eq. ©, © we find 



vWA^T)) = ^ MdX^{g2)), (166) 
^ Pn{dXi{l)) ^ g2^ Pn{dX'i(g2)) n=l,2,3,..., (167) 
a„(dA«(l)) = <?2a„(dA«(g2)) n = 1,2,3,.... (168) 



Hence by comparing Eq. (jl66p . (jl67p and (jl68p with the chain coefficients in Eq. (jl38p 
and Eq. (jl59p , we see that the chain mapping coefficients are invariant under any choice of 
52 > 0. □ 

Corollary 4 The generalised mapping Eq. (jl22p reduces to 

00 

H^Hs + ^/M^A{bo + bl) + Y, 52a„(0)6t fe„ + 52 V/3«+i(0)(6l+i6„ + /i.e.). (169) 
w/ien 9 = 0, yl2 = 0. We also find 

(170) 

(171) 
(172) 

This is just the known result of found in |19j and |18) . 

Proof Follows from setting q — 0, yl2 = in theorem ([T5|) and simplifying the resultant 
expressions. □ 



M"(x) 


J{x) 

1 


TT 


9'{x) 




max (0)} 


= {0,l^max/g2} 
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Given that the couphng of the chain elements is excitation number preserving, the elemen- 
tary excitations of the chain can be viwed as particles hopping on a Id lattice. We therefore 
make the following definition. 

Definition 17 We shall refer to the transformation described in corollary (|4]) as the particle 
mapping. 

Remark 4 lemma (fTUl) clarifies the role of the constant g in [TH] as the g2 constant defined 
in this article is equivalent to the g constant in their article for the particle mapping. 

Corollary 5 When q = 1, A2 = 0, the generalised mapping Eq. (jl22p reduces to 

H^Hs + ^/M^AXo + f2 LV/3n+i(l)^n^n+i + 92^^ + . (173) 

where Xn and Pn are position and momentum operators, Xn ■— {bli + 6„), P„ ■= i{bn ~ 
bn)/2 n = 0,1,2,.... We also find 

g\x) = (174) 

M^{x) = (175) 

TT 

(l)} = {0,^L./52}. (176) 

Proof Follows from setting q — 1, = in theorem (I16|) and simplifying the resultant 
expressions. □ 

Given that the coupling of the chain elements in Eq. (I173P resemble that of springs obey- 
ing hooks law, the elementary excitations arc phonons such as in solid state physics. We 
therefore make the following definition. 

Definition 18 We shall refer to the transformation described in corollary ^ as the phonon 
mapping. 

Remark 5 In light of definitions ([T7| and ([T5)) , we note that Eq. (jl22p interpolates between 
the two solutions. 



Definition 19 We call mth embedding to the new system- environment interaction pro- 
duced when the new system is composed of the initial system plus the first m sites of the 
chain formed by the environment in the chain representation. The new environment is 
formed by the remainig chain elements. See figure (c). 

As we will see in section 13.41 the chain coefficients converge for a wide range of spectral 
densities, and hence all the specific features of an environment appear in the first sites of the 
chain. Consequentially, these can be progressively (or directly, all in one go) absorbed into 
the system by making an embedding; to reduce the complexity of the effective environment. 
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Corollary 6 The generalised mapping in terms of Jacobi matrices is: 



H = Hs + ^/Mq) M^'-^'bl + e-'^-bo) 



(177) 



,-21^2 



where 



(jidX") ~ |l) bo + /i.e.] + hi (jidX^) + |l) b, 



(6„, 6„+i, 6„+2, &„+3, . . .)^, n = 0, 1, 2, ... , 



(178) 



and f is understood as the transpose plus h.c. of the elements of the vector. After the mth 
embedding the Hamiltonian is 



+ 



-2iA2 



h.. 



m = 1,2,3,... , 

where the Hamiltonian of the system part of the mth embedding, Hgq^ is given by 



m-2 



J2 VP^^) {q{e''^'bibi+i + e-2'-^^&„6„+i) + (fet 6„+i + b^bl^,) 

Tl = 

m— 1 ^ 2\ 

+ E |""('^) - J i^'^^'bibi + e-2^^^6„6„) + («„(g) + I) 6t 6„. 
For the particular cases of the particle and phonon mappings, this reduces to 



(179) 
(180) 



(181) 
(182) 

(183) 



H = Hso + ^/ l3rn{dX0){blbr,,-l + k.C.) + hljrn{dX°)hrn 1,2,3,. 



H = Hsi^ + ^l3rn{dX^)X^X^^i + -X^J™(dAi)X„ + -P^P™ 



(184) 
(185) 



m= 1,2,3,... , respectively. WhereXn := (X„,X„+i,X„+2, ■ • 0^; P" := (Pn, Pn+i, Pn+2, ■ 
71 = 0,1,2,... . 

Proof Eq. (I177p follows from writing Eq. (jl59p in terms of the Jacobi matrix of the measure 
c?A' and Eq. (jl79p then follows using the identity JidX*^) — Jm{dX'^) which is a direct 
consequence of corollary (jSj. □ 

Definition 20 We call the nth partial spectral density J„(a;) to the spectral density which 
describes the system- environment interaction of the nth embedding. We call the initial spec- 
tral density Jo{uj) such that Jo{uj) = J{uj). See figure (c). 
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(a) CwwvqJw 

w w 



(b) 



System Environment 



Om-1 




New System 



New Environment 



Fig. 1 (a) The initial system-environment Hamiltonian before tiie mapping has been performed: 
the system (depicted by the gray ball in the centre) couples directly to the degrees of freedom of the 
enviroment (orange balls), (b) The system-environment Hamilatonian after the chain mapping has 
been performed: the environment has been mapped onto a semi-infinite chain of nearest neighbour 
interactions where the system now only couples to the first element in the chain, (c) The system- 
environment Ifamiltonian after the mth environmental degree of freedom has been embedded: the 
new system, interacts with the new enviroment via the mth partial spectral density Jm (x) . 



3.2 Connection between the phonon mapping to previous work and the sequence of partial 
spectral densities 

Theorem 17 The sequence of partial spectral densities in |21j are generated by 

JrM = ^71^^^ "^-^ (186) 

(P„_i(dAi;a;2)i£(£^^ _ Q„_i(dAi; w2))2 + j2(^)p2_^(^^i.^2) 
n — 1, 2, 3, . . . , where g2 — I- 

Proof As observed by Leggett [3S] , the spectral density of an open quantum system can be 
easily obtained from it's propagator Lo{z). The authors of [21] have developed this to find 
a continued fraction representation for the case of the mapping of the propagator presented 
in their paper as follovifs. 

Lo{z)^-z^-wo{z), (187) 
Mz)^ ^^^^^2 ' (188) 

02 _ .2 fli 



2 _ ^2 
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where 

2 



Df^ = - / dujJniuj)u} n = 0, 1, 2, . . . , (189) 



2 

^22 =—-/ dt^J„(w)w3 n = 0,l,2,.... (190) 

Alternatively, we note that we can write the continued fraction Eq. pSSp as a recurrence 
relation 

^"(^^ ^ ^^^^-^^ITTTi) - = 0'i'2.---' (191) 

where u'^(z) := — i(j„(z). From 21;, we have an alternative expression for w„ in terms of 
the partial spectral densities though the relation. 



u;„(z) = i Tdo.:^^ n = 0,l,2,..., (192) 



TT 



which by a change of variables and taking into account the definition of uj'^ in Eq. (|191l) . 
can be written in the form 

w'^{^z) = - r du^^^^ n = 0,l,2,.... (193) 
Furthermore, via a change of variables the integrals Eq. (|189p and (|190p can be written as 
Dl = -l dujJniV^) n = 0,l,2,..., (194) 

(195) 



1 f°° 

^n+l = —^ dujJn{\^)uJ 71 = 0,1,2,.... 



We note that in [2T] the support of the spectral densities corresponds with the domain of 
the spectral densities defined in definition (jl6p and hence integrals Eq. (jl93|) . (jl94l) and 
(|195p are zero outside of the domain of Jn{x), n — 0,1,2, . . ., therefore we can change the 
upper infinity limit of the integrals by Now let us define the set of measures 

djnit) ^ dt%{t) 0,1,2,..., (196) 

7nW = ^^^ n = 0,l,2,.... (197) 

This definition of the measure has some important consequences: 

1) From Eq. (|195p we note that ^^'^ the first moments of the measures djn{t), 

C\{djn)^nl+i n = 0,1,2,.... (198) 

2) From Eq. (|193p and definition (O we see that w^(yi) is proportional to the Stieltjes 
transformations of the measure d-jnit) 

wUV^) = DlSUz) n = 0,1,2,.... (199) 
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3) From Eq. (jl94p we see that the zeroth moments of the measures d'yn{t) are unity 

C^o(rf7n) = l n- 0,1,2,.... (200) 
We are now able to re-write Eq. (jl9ip in the form 

5„+i(z)A2+i =z-Ci(d7„)--l^, n-0,1,2,..., (201) 

where we have used the short hand S'm^(z) —: Sn{z). By comparing this recursion relation 
with Eq. (fTO)) and definition ([5]), we deduce that D^^idjn+i is the secondary measure 
associated with (!"/„, for n = 0,1,2,... . We can also identify a sequence of normalised 
secondary measures. Noting that Ca{d'yo) = 1 from Eq. (|200p . definition © tells us that 
the sequence of secondary normalised measures starting from d'ye is 

d7o, Dldji/Co{Dldji), Dld^2/Co{Dld^2). ■ ■ . , Dldj^/Co{Dldj^), .... (202) 

However, 

Co(A'rf7n) - DlCo{d-fn) ^Dl n = 1, 2, 3 . . . , (203) 

so the sequence of normalised secondary measures is d'yo{t), d'yi{t), d'y2{t), d'y^lt),... . 
Taking into account lemma ([3]) and theorem ^ we see that 

Dl ^ d„^i ^ (3„{djo) n= 1,2,3.... (204) 

Due to corollary ([T|) and Eq. (|200p we can also write Dq in terms of /?o, 

Dl = /3o(d77o) (205) 

where 

drto DldjQ. (206) 

We can now construct a sequence of beta normalised measures from drjo, denoted by drjo, 
dr]i, dr]2, .... From Eq. (|205p and (|206|) we see that dijo satisfies Eq. (jlOip for ^70, hence 
lemma ^ tells us 

Vnix) ^ Pnid7^o)%{x) n = 0,l,2,.... (207) 

Taking into account that /3„((i?7o) = /3n('^7o) = 1,2,3,... due to lemma ([7]) and Eq. 
((206)) . from Eq. (|207)) and (ITWl) we gather 

dr]n{t) ^ dtfj„{t) n = 0,l,2,..., (208) 
fy„(t) = :M^ n = 0,1,2,..., (209) 

TT 

hence we note that 

d-nQ^dXK (210) 
Substituting this into Eq. (|M1) gives us Eq. (|186p . □ 

Theorem 18 T/ie mapping described in |21| zs a special case of the generalised mapping, 
theorem (|16p . namely it is equivalent to the phonon mapping. 
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Proof The mapping in 21 is 



H^Hs- DosX, + J2 ( -DnXnXn+i + + J ^ (211) 



n=l 



where Hs ^ £^ + V{s) + AV{s). See [21] for full details. From Eq. (IT95)) and (0) we see 
that we can write Sl^ as 



f2l+^ ^ aoidrjn) n = 0,l,2,.... (212) 



Hence, 



122^1 =a„(d7]o) =an(rfAi) =a„(l) n = 0,l,2,..., (213) 

where we have used Eq. (ITU5)) followed by Eq. (I^TU| and From Eq. ([TM|) . and dH) we 

see that 

Dl = Po{dVn) 0,1,2,.... (214) 

Hence using Eq. ([T09l) followed by Eq. ((2101 and (IT55|) . we find 

Dl = l3,,{d'no)^Pn{dX^)^Pn{l) 71 = 0,1,2,.... (215) 

If we let v4 = s in Eq. (|173l) and make the change of variable rt — > rt + 1 in the summation 
and relabel the position momentum variables such that (X„,P„) — > (X„_|_i,P„+i) n = 
0, 1, 2, . . . , from Eq. (PH)) and ((^15)) we find 



OO 



H^Hs + DosXi + i^DnXnXn+i + ^Xl + ip^^ j . (216) 

Now if we redefine the position and momentum variables by performing a phase rotation 
Xn -> X„(— 1)", P„ P,i(— 1)" to get a minus sign in the interaction terms and let 
iJs - ^ + y{s) + Z\V^(s) in Eq. (dZSl), we find Eq. (gill). □ 

Remark 6 We note that the three term recurrence relations Eq. ([S]) for measures with gaps 
(i.e. measures d^[x) — fi{x)dx where the support of fi.{x) is not one section of the real line, 
but rather has disjointed support regions i.e. / —[a, h] U [c, d\ where h < c), still hold. Hence 
the generalised mapping (and therefore the phonon mapping) is still valid. However, the 
chain mapping presented in [21] is not valid under these conditions because the Stieltjes 
transformation of a measure, definition ([T]) is not valid anymore and hence one cannot 
calculate the chain coefficients from the sequence of partial spectral densities. In this sense, 
the phonon mapping presented here is a more general result. This is an important difference 
because if a spectral density has a gap in its support, then the corresponding measure also 
has a gap. There are physical systems (such as photonic crystals and diatomic chains) which 
have these properties. 

Corollary 7 In the phonon mapping case, the sequence of partial spectral densities is given 
by 

JrM = TITT-^, '^-^ (217) 

n — 1, 2, 3, . . . , where 52 = 1- 

Proof Follows directly from theorems ([T7| and (|18p . Also see section 16.11 in the appendinx. 
□ 
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3.3 Sequence of partial spectral densities for the particle mapping case 

Theorem 19 In the particle mapping case, the sequence of partial spectral densities is given 
by 

Jn{^) = 77T<r-, (218) 

(P„_i(dA";c.)^^^-Q„_i(dA0;^))2 + j2(^)p2_^(dA0;^) 

n = 1, 2, 3, . . . , where 52 = 1- 

Proof The proof will start by defining a set of initially independent Hamiltonians. Each one 
of these will represent some system-environment interaction of a similar type discussed in 
this paper. When they are in this form, we can easily find an expression for the spectral 
density representing the system-environment interaction. We then put constrains on the 
Hamiltonians in such a way that they are no-longer independent from one another. Finally 
we show that these spectral densities correspond to the sequence of spectral densities for 
the particle mapping. 

Let us start by defining a set of independent Hamiltonians 

H,n — Hsm+ / 9{x)amxalnxdx + / hm{x){al,^xb bl,_i)dx (219) 



"'0 



m = 1,2,3 ... , where [ajnx,aiy] = SmnS{x-y) and [bm-i,bl_j^] = 6„in rn,n = 1, 2,3, . . . 

&„i-i, i>\n-i for the mth Hamiltonian Hm belong to the system which has Hamiltonian 
We note that the mth hamitonian H„n has a spectral density J™ (a;) given bjQ Eq. (ina 



J-^^) = ^hUg^\u)fj^. (220) 

Now let us set g{x) = g2X. Eq. (|220p gives us hf^^{x) — g2J"^{g2x)/Tr. We now define the set 
of measures 

dd,n{x) ^hl^{x)dx ^ g2J"'{92x)dx/n m = 1,2,3,..., (221) 
and the new set of creation and annihilation operators through the unitary transformation 

6t„„=/ dxU::'{x)al, n-0,1,2,..., m = 1,2,3,... , (222) 

Jo 

where U^{x) = hjn{x)Pn{ddm] x) n = 0,1,2,... ,m = 1,2,3, ... . Using the orthogo- 
nality conditions of orthonormal polynomials definition ([3]), we find the inverse relation to 
be 

00 

<x = Y.^n{x)bln m = 1,2,3,.... (223) 

n=0 



^ We have used a superindex here rather than a subindex to denote the spectral densities so as 
not to confuse them with partial spectral densities, as at this stage we cannot identify them as 
such. 
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Substituting this into the RHS of Eq. (j219l) and using the three term recurrence relations 
Eq. ([5]) and orthogonahty conditions Eq. (jS]), we find that we can write Hm as 

= Hs„^ + y/MdA^iblno^m-i + brnQbl_i) (224) 

oo 

+ h.C.) + an{ddm)b]nnbmn (225) 

n=0 

m = 1, 2, 3, . . . . We note that at this stage, the set of spectral densities {^™(a;)}m=i are 
independent and undefined. This freedom allows us to let the set of measures {d'dm]m=i 
be a sequence of beta normalised measures generated from the measure which from 
Eq. (|154l) and (jl75p we see is given in terms of another spectral density J{x). Now the 
spectral densities {J'"(a;)}^^i are fully determined by J{x) through the definition of a 
beta normalised sequence of measures, definition (|T3| . Hence using Eq. (l99l) and (j22ip we 
get 

^"(W) = 7TT^, — (226) 

(P„_i(dAO;c.)^£(^ - Q„„i(dA0;c.))2 + j2(c,)/'2_^(^^o.^) 

n= 1,2,3,... , where where 52 = 1- Now we have to show that J" is the nth partial spectral 
density for the particle mapping Hamiltonian Eq. (jl69p i.e. the interaction described by Eq. 
HHH). First let us define Hsm in Eq. (gM]) as Hsm = Hgo^ m = 1, 2, 3, . . . and note that 
using Eq. and (fTUg)) . we have 

/3o(d^m) = /3m(dA") m= 1,2,3,..., (227) 
a„(di9m) = a«+m(c?A") m = 1, 2, 3, . . . , n = 0, 1, 2, . . . , (228) 
/3„+i(dz9™) = /3„+i+™(dA") m= 1,2,3,... , n = 0, 1,2,... . (229) 

Hence Eq. (|224p now becomes 

ff„ = Hs« + v/An(rfA")(&]„o&,„_i + 6™o6L_i) (230) 



+ v//5n+i(dAO)(6L(„+i_„)6„(„_™) + /i.e.) (231) 

n—m 

+ a,„(dA°)6|„(„_„^)6™(„_™) (232) 
m — 1, 2, 3, . . . , or by writing it in terms of a Jacobi matrix, 

H,^ = Hso^ + V/3™(rfA")(6j„o6„,_i + h.c.) + bJ„o J'„.(dA°)b™o m = 1, 2, 3, . . . , (233) 

where b,^„(, := (fe^, fotii- ^L- ^Ls' • • •)> b™o := {bmo, bmi,b,n2, b,n3, ■ ■ -V m = 1,2,3,... . 
By comparing Eq. (|233l) with Eq. (|184p . we see that we have generated a chain mapping 
with identical coefficients, hence the spectral densities { J™(a;)}^^]^ are the partial spectral 
densities for the particle mapping. □ 
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3.4 Partial spectral densities sequence convergence 

Definition 21 We say that the chain mapping for some q and a particular spectral density 
J{x), will belong to the Szego class if the measure d\'^{x) ~ M''{x)dx satisfies 

^,.a.{q) \nM'}(x)dx 

^ ' > -00. (234) 



„_(g) \/ i^max (?) x)(x X^in^ (?)) 

Remark 7 Examples of spectral densities which for any q do not belong to the Szego class, 
are those which have gaps in their domain and those with unbounded support. 

Theorem 20 If for some q-chain mapping for spectral density J{uj) belongs to the Szego 
class, then the sequences ao(g), ai(g), 02(9), • . • and ^^{q) , I3i{q) , I32{<l) , ■ ■ ■ converge to: 

hm a„(9) = (235) 

n— >oo 2 

hm /3„(g) = ^ — '—. (236) 

n— >oo 10 

Proof Follows from shifting the support region of the Szego theorem in [TS] . For the original 
theorem see [3^. □ 

Corollary 8 If for some q-chain mapping for spectral density J{uj) belongs to the Szego 
class, the tail of the semi-infinite q-chain mapping tends to a translational invariant chain. 

Proof Follows from theorem (^01) and Eq. (|159p . □ 

Definition 22 The moment problem for a measure dfi is said to be determined, if it is 

uniquely determined by its moments. 

Theorem 21 If a measure dfj, has a finite support interval I, then its moment problem is 
determined. 

Proof See [17]. 

Theorem 22 If for some measure dfi{x) with finite support interval I the limits 

lim a„(d^) = (237) 

n~^oo Z 

hm Mdfi) - (238) 

exist, then the sequence of beta normalised and normalised secondary measures generated 
from dv and dfj, respectively, converge weakly to 



lim P„(.) ^ ^(-~-)('--\ (239) 



8^{x-a){b-x) 

hm ^in{x) = — jT . (240) 

n-i-oo 7r(o — a)^ 



Furthermore, this is the only converging limit which exists. 
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Proof First we will show that the limit exists and is unique by construction, then we will 
show that it corresponds to when the limits Eq. (j237l) and (j238p are accomplished. 
Substituting Eq. §0^ into Eq. ([Ml), we find 

^idp;x) = 2[x- Ci(d^)] - ^^^^J^y^'^^ (241) 

If the limit exists, then there must be at least one solution to p{x) — Ap,(x) where dp{x) — 
p{x)dx is the secondary measure associated with the measure dp{x) — p,{x)dx and A > Q. 
From definitions (ITT|l and ([7]) we see that f{dp] x) = A(p{dp; x) and hence from Eq. (|24ip 
we see that 

^idp;x)^^-^. (242) 
Thus substituting this into Eq. ([5^ . and solving for p{x) we find 

Taking into account the definition of a measure, definition ([T]), we see that if the limit exits, 
then it must be bounded. Moreover, it must belong to the interval centered at C'i{dp) and 
of length 4^/A. Making the change of variable x — Ci{dp) = 2\J~At gives us 



2J\ - t2 

dpit) = p{t)dt = — dt, (244) 

TT 

with support interval [-1,1]. We note that at every step of the proof we have used the most 
general conditions and found the most general solution, hence if this solution exists, it must 
be unique. We now can check that Eq. (j244l) exists by direct substitution. Using definition 
([7]) we find that the Stieltjes Transform of Eq. (|244|) is S^{z) ^ 2[z - \J z"^ - 1] and hence 
from definition (|lll) we find ^(dp\ x) — Ax. Thus using Eq. ([M]) we get p{x) = /i(a;)/4, hence 
the limit exists. By performing the change of variable t = (2/(6 — a))y + (6 + a)/(a — b), 
we shift the support of Eq. (|244l) to the general case [a, b] and the measure is now given by 



Eq. (12401) . 

Now we will proceed to show that if Eq. (j237p and (j238p are satisfied, then the sequence of 
normalised secondary measures converge weakly to Eq. (j240p . By taking the nth moment 
of the measures in Eq. ([^^ and taking into account lemma (1^^ for n = , we have 

C'nidpn+l) = [C2{dpn) ^ Ci{dpn)^] C„(d//„+i) n, TO = 0, 1, 2, . . . . (245) 

Writing Eq. (|24l) for a sequence of measures followed by substituting in our expression for 
Cnidpn+i) using Eq. (|245p . we find 

n-l 

{4 - (4?) = - - 5] (c^ - {cAf) d^'c^., n, s = 0, 1, 2, . . . , (246) 

where := Cn{dps) n,s — 0,1,2,.... Let us define the limit /„ = lim^-^ooC^- We can 
now draw the following conclusions 

1) Due to corollary ([5]), we have Co(rf/i„) = 1 n = 0,l,2,.... Therefore = 1- 
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2) From Eq. ([57]) and © we see that 

a„(dMo) = Md^^n) = = c?. (247) 

Hence taking into account assumption Eq. (|237p . we have li = {a + b)/2. 

3) Noting the definition of d„ in theorem theorem (|S]) tehs us 

/3„+i(dAio) = 4'-(c?)' n = 0,1,2,.... (248) 

Hence taking into account assumptions Eq. (|237p and (j238p . and Eq. (|247p we have 
^2 = (5a2 + 6a6+562)/16. 

We note that for any s, all the moments in Eq. (|246l) are fully determined by the starting 
values Cq, cf , and c|, hence we conclude from the above points 1), 2) and 3), that under the 
assumptions Eq. (|237p and (|238p . all ^„ 7i = 0, 1, 2, . . .are finite and determined by 

n-l 

{h - ll) In = ln+2 - hln+l - {h - if) J] ¥n-J n, S = 0, 1 , 2, . . . , (249) 

with starting values Iq — 1, li — {a + b)/2, and I2 ~ (5a^ + 6ab + 5b^)/16. 
For the case of the normalised measure Eq. (|240p . we conclude that — cf^^ n,s — 
0,1,2,... since its secondary normalised measure is equal to itself. Hence by denoting m„ — 
n, s = 0, 1, 2, . . . we can write Eq. (|246p for this measure as 

n-l 

(m2 - (toi)^) to„ = m„+2 - rriimn+i - (rn2 - (™i)^) mjinrin-j (250) 

n,s — 0,1,2,... . By direct calculation of the moments of Eq. (I240p . we find that mg = 
^0, rni — /i, m2 — h, and hence by comparing Eq. (j249p with Eq. (I250p we conclude 
that m„ — In n — 0,1,2, ... . Thus using lemma ([2T|) we conclude Eq. (|240l) under the 
assumptions Eq. (^771) and ((255)) . For Eq. (USD), we note that Eq. dTUIH) and Eq. ([TUS]) tell 
us lim„_j.oo j^n(a;) = limn-i-oo /3n('^Mo)A*n(a^)- Hence from Eq. (I238P and (|240p we conclude 
Eq. ([^ . □ 

Definition 23 W^e wiH call terminal spectral density Jt^i^) the spectral density to which 
a sequence of partial spectral densities converge to weakly if such a limit exists: Jt(w) = 
lim„_>oo Jn(w). 

We are now ready to state our fourth main theorem: 

Theorem 23 // for the particle or phonon mapping the spectral density J{lu) belongs to 
the Szego class, the sequence of partial spectral densities converge to the Wigner semicircle 
distribution |37| 

and the Rubin model spectral density jlj 

^ (252) 

respectively. 
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Proof From Eq. (^0^1) and theorem ^ we gather that if Eq. ((^57)) and are satisfied, 

then for the phonon case if we choose 52 = 1 

v/(c.^-a)(b-c.2) 
Jt(w) = . (253) 

From Eq. (jl76l) we have that a — and b — Now taking into account theorem (I^Ul) . 

we get Eq. (j252l) . Proceeding in a similar manner, we find Eq. (I25ip . □ 



4 Examples 

4-.1 power law spectral densities with finite support 

The the widely studied power law spectral densities is pHlfT] 

J{x) = 27701^1-" x", (254) 

with domain [0,aJc] and s > —I. Let us start by calculating the sequence of partial spectral 
densities for the case of the particle mapping. 
From Eq. (ITTO)) and (PMl) we have 

M°(x) = 2aujl-~'x''. (255) 

For simplicity, we will scale out the uJc dependency and show how to put it back again 
afterwards. Let us start by defining the weight function mQ{x) := u!cM'^{xu!c) = 2au!Qx'^ 
with support interval [0,1]. From lemma ((H]), we find that rn'^{x) = M^{xujc)/'^c n — 

1, 2,3, where m^^{x) and M^{x) are the sequence of beta normalised measures generated 

from mQ(x) and Af°(x) respectively. Now let us define the weight function 

m[J(a;) = x" (256) 

with support interval [0,1]. Given that this new measure is proportional to mQ(x), from 
lemma (1131) we conclude that it's sequence of beta normalised measures 'm'^{x) n = 
1,2,3,... , are equal. Hence we have M^{x) = u!crn'^{x/ujc) n = 1,2,3,... . Thus taking 
into account Eq. (|170p and (|221[) . setting 92 1 we conclude 

J„(cj) = 7rM°(w) = Wc7rm° (ljK) n = 1, 2, 3, . . . . (257) 

The real polynomials orthogonal to the weight function mQ(a;) are Pn{x) ■— Pn'^\2x — 
l)^/n + s + 1 n ~ 0,1,2,... ; which are normalised shifted counterparts of the Jacobi 
polynomials Pn°''^\x) n = 0, 1, 2, . . . . The reducer for the case s > is given by theorem 

r-1 



ip{dm^;t) 



In ( ]+s I x'^^ In 



t — X 



t 



dx 



(258) 
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which has analytic solutions when a particular value of s is specified. For example, the first 
three in the sequence for the ohmic case (s = 1) are 

= 2(.2.^ + [l + .ln(i^F)' 



^2{x) = . . „ Z _ . ,._^..2 ^ (260) 

rhlix) 



47r2(2- 3a;)2a;2 + [l - 6a; + (4 - 6x)xln(i^)] ^ ' 

6x 



367r2a;2(3 - 12a; + 10a;2)2 + [30a; - 16 + (18 - 72a; + 60a;2)(l + a;ln(i^))] 
The chain coefficients are calculated in [l8l to be 



""("^(-)-yO+ (,s + 2n)(2V. + 2.) ' - = M,2,..., (261) 



V/3n+i(0)(.) - ^ 2 + 2n)(3 + ,s + 2n) V 1 + + 2n "^O'^'^,..., (262) 

(263) 

where the (s) is to remind us of their s dependency. System environment coupling coefficient 
is dH] 



Vm=^c^j^- (264) 

Now we will find the sequence of partial spectral densities for the case of the phonon 
mapping. From Eq. (|175p and (|254p we have 

M\x) = 2aujl-'x'/^, (265) 

with support [0,a;2]. Proceeding as in the particle mapping case and taking into account 
Eq. (pngi) . we find 

J„(w) = 7rAfi(a;2) =tj27rmi(tj2/^2) n 1, 2, 3, . . . , (266) 

where rhl^{x) n — 1,2,3,..., are the sequence of beta normalised measures generated 
from 

mJ(x)=x''/^ (267) 

with support [0,1]. Let us denote 'itiq{x) and ttiq^x) by mjj,(a;) and ifiQ^^x) respectively 
to remind us of their s dependency. By comparing eq. (j267p with eq. (|256p , we have that 
"^Os(^) — ™os/2(-^) f*-"^' all s > — 1 and hence from Eq. (12571) and (|266p we see that for the 
Spin-Boson models, there is a simple relationship between the partial spectral densities of 
the particle and phonon mappings for different s values. For the same example as in the 
particle case (s — 1), we need to evaluate m° (a;) n = 1, 2, 3, . . . , for s = 1/2. The first one 
in the sequence is 

m](x) = — . ^, (268) 

' 3(7r2a; + (2-2Vitanh~i(Vi))2)' ^ ' 
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where tanh~ {x) is the inverse hyperbohc tangent function. We can readily calculate the 
chain coefficients from the particle example. By comparing the expression for the a„ and 
/3„ coefficients for the weight functions for the particle and phonon mappings, we find 



an{l){s) =uj,an{0){s/2) n = 0,1,2,..., (269) 

= Lo,Vl3n+ims/2) n = 0,l,2,..., (270) 

and system environment coupling term to be 

\/^=2c.,y^. (271) 



For both the particle and phonon mappings, it is easy to verify that the chain coefficients 
and sequence of partial spectral densities will converge because Eq. (j234p is satisfied in 
both cases as long as s < oo. We also see that the sequence of partial spectral densities 
calculated in the above examples converge very rapidly to this limit after about the 3rd 
partial spectral density. 

4-2 The power law spectral densities with exponential cut off 
The power law spectral densities with exponential cut off is pHlfT] 

J{x) = 2T:aujl-'x'e-'-'/'^% (272) 

with domain [0, oo) and s > — 1. Let us start by calculating the sequence of partial spectral 
densities for the case of the particle mapping. 

From Eq. ^70)) and (j??^ we have 

M^{x) = 2aujl-'x' (273) 

Let us define the measure m'^{x) :— uJcM'^{xuJc) — 2aujQx'^e^^ with support interval [0, oo). 
From lemma (fT2)) . we find that m^^{x) = M^{xuJc)/uJc n = 1,2,3, .. .. where rn'^{x) and 
M°{x) are the sequence of beta normalised measures generated from m^{x) and M'^{x) 
respectively. We will now define a 3rd measure by m°(a;) = x^e~^. We note that it is 
proportional to the measure rrfiix) and hence lemma (jl3p tells us that its sequence of beta 
normalised measures are equal. Thus we have the relation 

Ml{x)=uj,ml{x/u:,) n = 1,2,3,.... (274) 

The real polynomials orthogonal to the weight function fh^{x) are called the associated 
Laguerre polynomials L'^{x) n — 0,1,2, . . . . Their normalised counterparts are Pn{x) '■= 
L!^{x) n\/r{n + s + 1) n = 0, 1, 2, . . . . The reducer in this case is given by theorem ((25)) 

(p{dm°;x) ^2 {s-t)t'-^e-Hn 
Jo 



t — X 



dt. 



(275) 
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which has analytic values when a particular value of s is specified. For example, if s is 
integer, we get 



k=s-l 



""Ex^x)- (s-fc-1) 



\x^ 



k=0 



(276) 



where Ejx{x) is the exponential integral function|38). 
From Eq. ([^^ followed by Eq. (E7il) we have 

Jn{u)^TiM^S^)^Lo^-Kml{uj/Lo^) n = 1,2,3,..., (277) 

and Eq. ([M]) tells us 

m°(x) = - ^^^4l n= 1,2,3,.... (278) 

For example, the first two in the sequence for the ohmic case (s = 1) are 
m^i^x) = , (279) 

g2a: _|_ j^2rj.2 _ 2xEx{x)e^ + x'^Ex{x)'^ ' 

— / \ , . 

" e2^(l - xf + n^x^ix - 2)2 - 2a;(2 - 3x + x^)Ex{x)e^ + x^{x - '2.)^Ee{xY ' ^^^^^ 

We also have analytic expressions for the chain coefficients. From [18], we have that the 
chain coefficients are 



V/3„+i(0) - i^c^{n + l){n + s + 1) n = 0,l,2,..., (281) 
q;„(0) = Wc(2n+l + s) n = 0,1,2,..., (282) 

with the system-environment coupling coefficient given by 



x//3o(0) =t^cv/2ar(s + l). (283) 

Similarly, we can calculate the partial spectral densities and chain coefficients for the phonon 
mapping case. 

Because the support interval is infinite, the Spin-Boson model with exponential cut off does 
not belong to the Szego class for either the particle mapping nor the phonon mapping as 
can be easily verified from Eq. (|234l) . This is reflected in the example above as the sequence 
of partial spectral densities do not convergeH. 

Remark 8 We note that for spectral densities for which the corresponding orthogonal poly- 
nomials are unknown analytically, we can easily calculate their coefficients using very stable 
numerical algorithms. See |19) or |18) and references herein for more details. 



* However, we do note that the ratio a„ {q) / (3n+i{<l) for g = and q = 1 does converge. This 
suggests that there is a universal asymptotic expansion for large n for the nth partial spectral 
density. 
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5 Summary and conclusion 

By developing the method of [T51IT5] , we have established a general formalism for mapping 
an open quantum system of arbitrary spectral density of the form Eq. (|119l) onto chain rep- 
resentations. The different versions of chain mappings are generated by choosing particular 
values of 4 real constants. This has also provided a very general connection between the 
theory of open quantum systems and Jacobi operator theory as the semi-infinite chains can 
be written in terms of Jacobi matrices as can be seen in corollary There has been a 
wealth of research into the properties of Jacobi operator theory [5Sj, and hence this opens 
up the theorems of this field to the possibility of being applied to the theory of open quan- 
tum systems. Likewise, the new theorem regarding Jacobi matrices (jl5|) developed here, 
could turn out to be useful in the field of Jacobi operator theory. 

There were two previously known exact chain mappings; the one of |18| which is the same 
as the particle mapping defined here (definition [T7|. and the mapping by [5T]. We point out 
that the phonon mapping derived in this article (jl8p has a wider range of validity (remark 
[6]) than the mapping of |21j and prove that they are formally equivalent in the range of 
validity of [H]; see theorem (IT5|) . 

The concept of embedding degrees of freedom into the system has been around for some 
time [Mini]: however, we see that in chain representations, there is a natural way of shifting 
the system-environment boundary, that is to say, there is a natural and systematic way of 
embedding degrees of freedom into the system one by one (or all in one go). To solve quan- 
titatively this problem, we have to first embark on finding the solution to an old problem in 
mathematics; an analytical solution to the sequence of secondary measures in terms of the 
initial measure, it's associated orthogonal polynomials and reducer; see theorem p^ . Not 
only does this provide a means to find analytical expressions for the sequence of spectral 
densities corresponding to the new system-environment interaction after embedding envi- 
ronmental degrees of freedom into the system for the particle and phonon mappings in the 
gapless spectral densities case (theorem and corollary ^ respectively), but it provides 
physical meaning to this abstract mathematical construct. 

Using convergence theorems of Szego and deriving the fixed point in the sequence of sec- 
ondary measures; we have combined these results to obtain a convergence theorem of the 
sequence of partial spectral densities for the particle and phonon cases, theorem (1231) (or 
equivalently, the sequence of secondary measures, theorem (HH)). What is more, because 
the criterion for convergence (definition \Tl\i is solely in terms of the initial spectral density 
once the desired chain mapping is chosen, the convergence criterion is readily applicable 
to a particular problem. Furthermore, we see that any unbounded spectral density will not 
satisfy this criterion. This is refiected in the sequence generated in the examples section 
for the case of the family of power-law spectral densities with exponential cut off as the 
sequence does not converge, section 

We give two examples where we can find explicit analytical expressions for the chain co- 
efficients and the sequence of partial spectral densities. These are the family of spectral 
densities used in the Spin-Boson model, which have spectral densities of the form a;* with 
finite support / — [0,uJc] and x^e^l'^" with semi-infinite support / = [0, -foo), where s > —1. 
Furthermore, we show how the partial spectral densities of the mapping for both phonon 
and particle cases are related for different families of the mapping as demonstrated by the 
identity mjs(x) = m°,^{x) and Eq. ((257)) and ([^7^ . 
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We note that when this convergence of the embeddings is achieved (i.e. inequahty (j234p is 
satisfied) , the part of the chain corresponding to the new environment has a very universal 
property: ah its couphngs and frequencies become constant. This is to say, they are trans- 
lationaUy invariant. This suggests a universal way of simulating the environment as all the 
characteristic features of the environment are now embedded into the system, as discussed 
for the mapping of [IH] in [22]. What is not so clear however, is what is the most effective 
chain mapping for simulating the environment. The general formalism of chain mapping for 
open quantum systems developed here, paves the way to answering these questions. 



6 Appendix 

6.1 Alternative proof for the sequence of partial spectral densities for the phonon mapping 

We note that by applying the same line of reasoning of the proof for theorem to the 
phonon mapping Hamiltonian Eq. (I173p . we can easily provide an alternative proof for 
corollary ([7]). This has the advantage of being an independent derivation of the results of 
|21j , but has the downside of not illustrating the connections between this paper and their 
results. 



6.2 Methods for calculating the reducer 

Theorem 24 If for a measure d^{x) — fl(x)dx, p-{x) satisfies a Lipschitz condition over its 
support interval [a,b], then 

.fidfi; x) = 2fi{x) In ( - 2 /' ^^^^l^J^dt. (284) 
\o-xJ 7„ t-x 

Proof From definitions ([TT|) and ([7]) we get 

^(dM;.)^lim2r ;--ffft - (285) 

Writing this as 

N ^ /"^ {x~t)fi{x) , „ f'' (x ~ t){p(t) ~ flix)) , 

^id,, X) . Inn 2 1 ^^-^dt + hm 2 ^J\% ^t (286) 

= hm 2 f l^-Z^dt - 2 r m^,t (287) 
2e2 \ /2(i) - /2(x 

e^Q^Ja \{x-tY+(^) X-t 

1) The first integral in line (|287p reduces to 

lim 2 t (^-^y(^) = lijn fL{x) \\n{{a - xf + e^) - ln((6 - xf + e^)] (289) 

= 2^(a;)ln ( ) . (290) 



li- / ( 7^^^ ) ^-'^^^r^dt. (288) 
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interval [a, b], the absolute value of the expression on line (I288P reduces to 



2) Now imposing the Lipschitz condition on ^i{x): \fi{x) — fi{t)\ < K\t — x\ for some K over 

(291) 

< lim / - — -dt (292) 



lim 



2e2 



(x - i)2 + 



\t~x\ 



dt 



b — x 



lim 2Ke 



arctan 



arctan 



= 0. 



(293) 
(294) 



Hence expression on line (|288l) vanishes. 



Lemma 11 // f{x) possesses a bounded continuous first derivative on its domain, then it 
also satisfies a Lipschitz condition on its domain. 

Proof See 

Theorem 25 If for a measure dfi{x) — fl{x)dx, p-{x) possesses a bounded continuous first 
derivative on its support interval I and p,{x) can be evaluated on the limits a, b, then 



ip{dfi; x) = 2 



fj,{a) In 



fiib) In 



fi'{t) In 



t — X 



dt 



(295) 



where fJ-'{t) is the first derivative of ^{t). 

Proof Follows from integrating Eq. (|284l) by parts and applying lemma (fTTj) . □ 



6. 3 Scaling properties of the sequence of beta normalised measures 

Lemma 12 Suppose we have two measures dv^(x) = i>^{x)dx and dv^ix) = iP{x)dx with 
support intervals bounded by Xa,Xb and a,b respectively. If they are related by 



V\x) 



v'^jx/X) 
X 



A > 0, 



(296) 



then their sequence of beta normalised measures have the relation 

i?^{x) = Xvl{x/X) 1,2,3,.... (297) 

Proof Using the definition of inner product, definition and relation eq. (|296p . we find 



{f,9)u^ = / fiXt)giXt)diy'{t). 



Hence we conclude 



Pnidv^;x)^Pnidv^;x/X) 0,1,2,. 



(298) 
(299) 
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Using this relation and definition ([B]), we find 



Qn{dv \ x) = 



n = 0,l,2,... . 



(300) 



Similarly, we have from definition PT|) 



ip{dv^; x) 



(301) 



A 



Finally, substituting Eq. (PMl) . pHI)) . and ((5UT|) into definition (US]) we arrive at Eq. (^^71) . 



Lemma 13 The sequence of beta normalised measures generated from a measure dv{x) — 
i?{x)dx are invariant under the mapping di'{x) — > \di'{x), A > while keeping the support 
interval I unchanged. 

Proof Follows with a similar line of reasoning as for lemma ((HI). □ 
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